Evidence for Ground- and Excited-State Efimov Trimers in a Three-State Fermi Gas 
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We observe enhanced three-body recombination in an ultracold three-component ®Li Fermi gas 
with large but unequal scattering lengths attributable to an excited Efimov trimer state near the 
three-atom scattering threshold. We find excellent agreement between the measured three-body 
recombination rate and the recombination rate calculated in the zero-range approximation where 
the only free parameters are the Efimov parameters and 77, . The value of k, determined by the 
location of the Efimov resonance we observe at 895 G also predicts the locations of loss features 
previously observed near 130 and 500 G [l|, Ql suggesting that all three features are associated with 
universal Efimov trimer states. We also report on the first realization of a quantum degenerate 
three-state Fermi gas with approximate SU(3) symmetry. 

PACS numbers: 21.45.-v, 36.40.-c, 34.50.-s, 67.85.Lm, 67.85.-d 



A landmark result of few-body physics is Efimov's solu- 
tion to the quantum mechanical three-body problem for 
low-energy, identicalparticles with a large s-wave scat- 
tering length a 0, 0]. Being independent of the par- 
ticular character of the two-body interaction. Efimov's 
predictions are universal. At resonance (a ±00), an 
infinite series of arbitrarily shallow three-body bound 
states (Efimov trimers) exist with a geometric spectrum: 
En = -(e^V^o)-" h^i^Jm (for n = 0,1,2,...). Here, 
So ~ 1.00624 is a universal constant, m is the particle's 
mass, and the entire spectrum is fixed by a single three- 
body parameter which is determined by short-range 
interactions [^. For large but finite a, the spectrum of 
Efimov trimer states are universal functions of a and . 
Further, all low-energy scattering obscrvables display a 
log-periodic dependence on a and the collision energy Q . 
These startling predictions ignited a search for Efimov 
physics in a range of physical systems 0, 0] • 

The first experimental evidence for Efimov trimers was 
obtained by observing both a resonant enhancement and 
interference minimum of three-body recombination in an 
ultracold gas of bosonic Cs near a Feshbach resonance Q . 
The Efimov effect in bosonic systems near an isolated Fes- 
hbach resonance is now becoming well understood follow- 
ing the observation of atom-dimer resonances Q , Efimov 
resonances between bosons with different masses [l^l and 



tests of universal log-periodic scaling |lll. Il2| 



The Efimov effect is also expected to occur for three 
fermions in distinguishable states (|1), |2), and |3)) if 
at least two of the three pairwise scattering lengths a.y 
(ai2, a23, and 013) are larger than the range of interaction 
(which, for atoms, is the van der Waals length .^vdw — 
\/ mCg / ft^ ) [El- For a Fermi gas of ^Li atoms in the three 
lowest spin states, investigated here, all three scattering 
lengths can be resonantly enhanced by three overlapping 
Feshbach resonances (see Fig [1] (a)) [l3| allowing for a 
study of the Efimov effect in a Fermi system which is 
deeply in the universal regime jcy | ^ ^vdw- For this ^Li 
system, the scattering lengths are typically unequal, and 
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FIG. 1: (color online) (a) Three overlapping s-wave Feshbach 
resonances in ^Li for states |1), |2), and |3) [l^. At high 
field, aij — > — 2140ao. (b) The binding energies (i?i2, E23, 
En) of the universal dimer states associated with the Fesh- 
bach resonances. The dashed lines (n and n') depict the bind- 
ing energies of the Efimov trimer states (not yet accurately 
known). Resonant three-body recombination occurs (arrows) 
when an Efimov trimer intersects the three-atom scattering 
threshold. The grey-shaded areas identify non-universal re- 
gions where E < E^dw = f? jmS^^Y^ or \aij\ < 2£vdw- (For 
®Li, ^vdw ~ 62.5 ao.) 

any number of them can be positive. This gives rise to 
an Efimov spectrum with a rich structure which is only 
now becoming understood theoretically IJ, |l5[ . Entirely 
new phenomena related to the Efimov effect have been 
predicted in this case (e.g. interference minima in the 
rate of atom-dimer exchange reactions) [l^. 

Beyond few-body physics, superfluid phenomena in de- 
generate three-component Fermi gases have recently been 
a subject of intense theoretical interest [l^i [13, [mH 
[20j . For fermions with ai2 ~ 023 ~ 013 < 0, there is a 
competition between trimer formation and several types 
of Cooper pairing [l6j , breached pair and unbreached pair 
phases [l8|, as well as a unique interplay between mag- 
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nctization and superfluidity [l^j 0] • Degenerate three- 
component Fermi gases open a window of opportunity 
to determine which pairing symmetries and phases are 
ultimately favored by nature as parameters are tuned. 

In this Letter, we report resonantly enhanced three- 
body recombination in a three-state Fermi gas with large, 
negative and unequal scattering lengths that can be ex- 
plained by an excited Efimov trimer state near the three- 
atom scattering threshold. We measure the three-body 
recombination rate for fields between 842 and 1500 G 
which is deeply in the universal region at high field (see 
Fig. mb)). The data is well described by a calcula- 
tion of the three-body recombination rate at threshold 
in the zero-range approximation for which the only in- 
put parameters are the known scattering lengths ai2, 023 
and ai3 [l3l | and the complex-valued Efimov parameter 
K* exp[z77*/so] which we determine. Here, 77* has been 
introduced to account for the nonzero width of Efimov 
states due to decay to deeply bound dimer states. We 
find that the value of k* predicts that Efimov trimers also 
cross threshold near 130 and 500 G. This strongly sup- 
ports the interpretation that the loss features previously 
observed [H, Q are related to the Efimov effect. Finally, 
we produce the first degenerate three-state Fermi gas in 
the SU(3) symmetric regime where 012 « 023 ~ 013. 

The first experiments with degenerate three- 
component Fermi gases investigated their stability 
against decay via three-body recombination [l|, [l] . These 
experiments studied a ^Li gas in the three lowest spin 
states (|1), 1 2), and |3)) and measured the three-body 
recombination rate as the scattering lengths (ai2, 023 
and 013) were tuned in a magnetic field. Enhanced 
three-body recombination was observed near fields of 130 
and 500 G where |ai2|, |a23| ~ ^vdw and |ai3| ~ 10i!vdw- 
Several authors used Efimov's theoretical framework 
to explain the magnitude and variation of the mea- 
sured three-body recombination rate between 100 and 
500 G despite its questionable applicability near the 
boundaries of this field range where |ai2|,|a23| ~ ^vdw 
14 1 21, 22, [2^. In this interpretation, the resonant 



enhancement near 130 and 500 G is due to Efimov-like 
trimers crossing the three-atom threshold near/in the 
non-universal regime. 

This same three-state ^Li mixture can be used to test 
universal predictions of Efimov physics with large scat- 
tering lengths since (1) |ai2|, |a23| and |ai3| > €vdw for 
fields above 660 G due to three overlapping Feshbach res- 
onances [l3l and (2) there is negligible two-body loss 0. 
The three-body recombination rate in the vicinity of the 
Feshbach resonances was measured in Ref. Q. However, 
the authors noted that the highest recombination rates 
were unitarity limited and therefore could not provide 
tests of universal predictions [24 1. 

Efimov's theoretical framework makes universal pre- 
dictions for three-body scattering observables in the 
threshold regime where the collision energy is the small- 



est energy in the system. For identical bosons with 
large negative scattering lengths, the characteristic en- 
ergy scale is set by the height of a barrier in the adiabatic 
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0.l58h^ /ma^ |2J]. In the 



three-body potential. 

three-state ^Li gas we study, where 012 < 023 < 013 < 0, 
we expect that for fields B > 875G (B > 960 G) 
which corresponds to scattering lengths |ai2| < 12, 250ao 
(|ai2| < 5, OOOao), a temperature T < 30nK (T < ISOnK) 
is required in order to compare our measurements to the 
calculated recombination rate at threshold. 

In order to reach such low collision energies, we adi- 
abatically release a two-state mixture of ^Li atoms pre- 
pared in a crossed optical dipole trap (as described in 
Ref. 0) into a larger volume hybrid magnetic/optical 
trap with smaller oscillation frequencies. We use two dif- 
ferent trapping geometries for the large volume trap to 
study atoms at T ~ 30 nK and T ~ 180 nK respectively. 
For trap A (B) , the combination of magnetic and optical 
(from a 1064 nm dipole trap) forces produce a trap with 
frequencies = 15(2) Hz, Vy = 0.242/BHz ± 1%, and 
ly^ = 12(1) Hz (z/^ = 33v^l + 1.4 x 10-3(5 - 842 G) Hz± 
3%, Vy = 2ly/l + 3.6 X 10-3(5 - 842 G) Hz ± 3%, and 
= 94(2) Hz) [2^. Here, the dependence of the fre- 
quencies on the bias field B (in Gauss) for trap A(B) 
arises from the nonzero field curvature of a large (small) 
set of electromagnets used to produce the bias field B. 

Starting from an incoherent 50/50 mixture of atoms in 
states |1) and |2), we create a three-state mixture with 
equal population in states |1), |2) and |3) by simulta- 
neously applying radio-frequency (RF) magnetic fields 
which drive |1) — |2) and |2) — |3) transitions at the field 
B. Each of the RF fields has a resonant Rabi frequency 
57/27r ~ 10 kHz, is broadened to a width of 10 kHz and 
applied for a variable time (chosen to be at least twice 
the observed decoherence time) between 10 and 100 ms 
depending on the bias field. 

We have confirmed (for the fields studied here) that 
the decay of each of the possible two-state mixtures of 
the three states is consistent with one-body loss due to 
background gas collisions. One-body loss for these ex- 
periments gives a 1/e lifetime l/F = 2.8 s. In the three- 
state mixture, three-body recombination involving one 
atom from each spin state can occur in addition to one- 
body loss. In this case, the number of atoms in each of 
the equally populated spin states N{t) evolves according 
to N/N = — F — L^{n'^) where L3 is the recombination 
rate constant and (n^) is the average value of the squared 
density per spin state. For a thermal distribution at tem- 
perature T in a harmonic trap. N{t) evolves according to 
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where v = (vxVyVz) ■ Empirically, the temperature of 
the gas remains approximately constant. This is consis- 
tent with the fact that the energetic atom and molecule 
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FIG. 2: (color online) The measured three-body recombina- 
tion rate constant L3. Resonant three-body recombination 
near 900 G is caused by an excited Efimov trimer crossing the 
three-atom scattering threshold (see text). The inset shows 
typical data for N{t) and T{t) which are fit to extract L3. 
In the inset the solid (dashed) line shows the fit to a model 
which does not include (includes) anti-evaporation. 

produced in a recombination event have a mean free path 
much larger than the sample size and exit the cloud with- 
out depositing energy. Our measurements are consistent 
with a model that includes a rise in temperature due 
to "anti-evaporation" 0] (see inset to Fig. [2]). However, 
lower values of are obtained if the temperature is sim- 
ply assumed to remain constant. 

To determine the recombination rate constant L^, we 
measure for each field B the number and temperature of 
the trapped atoms as a function of time by in situ ab- 
sorption imaging. Wc show an example of data recorded 
at a particular field value in the inset to Fig. [51 The 
number evolution at each field is fit with an analytic so- 
lution to Eqn.[l] using L3, the initial number (iVo) and 
temperature (T) as free parameters. The recombination 
rate constant for fields between 834 and 1500 G is shown 
in Fig.[2]for data sets taken with traps A and B. The tri- 
angles (circles) correspond to trap A(B) and give L3 for 
a cloud at a temperature < 30 nK (< 180 nK) and an ini- 
tial peak density per spin state no — 5 x 10^ atoms/cm'^ 
(no ~ 5 X 10^° atoms/cm ). The error bars give the sta- 
tistical error from the fit and the uncertainty in the trap 
frequencies added in quadrature. A systematic uncer- 
tainty of ±60%, which arises from our uncertainty in the 
absolute atom number (±30%), is not included in these 
error bars. As wc vary the field, the scattering lengths 
arc tuned by broad Fcshbach resonances at 690, 811, and 
834 G and the observed value of L3 changes by several 
orders of magnitude. In addition to a smooth variation of 
L3, resonant three-body recombination is observed near 
900 G in both data sets. 

In order to compare our data to zcro-tcmperaturc cal- 



culations of the recombination rate, we require that the 
temperature be in the threshold regime and that the 
data is unaffected by the unitarity limit As argued 
above, the threshold regime should be obtained for the 
30 nK (180 nK) gas for fields above 875 G(960 G). For 
a gas at a low but nonzero temperature T, the recom- 
bination rate constant for three distinguishable fermions 
with equal mass is unitarity limited to a maximum value 
-/>3max = a/ikBTf whcrc a = ^^m^T^ h"^ / 0, 111- 
For a thermal gas with T = 30nK (T = 180 nK), 
^3max = 8 X lO'^^ cmVs (Lgmax = 2 X IQ-i^ cmVs). The 
measured values of L3 are < (1/I0)i3,nax for the 30 nK 
(180 nK) data for fields above 907 G (975 G). These sub- 
sets of the data should be in excellent agreement with 
zero-temperature calculations of the recombination rate. 

Recently, Braaten and co-workers calculated the three- 
body recombination rate constant in the zero-range ap- 
proximation for three fermions in distinguishable spin 
states at threshold [l4|. By numerically solving a gen- 



eralization of the Skorniakov— Ter— Martirosian (STM) 
equation, they calculate the recombination rate constant 
for the case 012, 023, 013 < 0. The only inputs to this 
calculation are 012, 023 and 013 (known for ^Li [l^l) and 
the two three-body parameters, and 77*, which respec- 
tively fix the spectrum and linewidth of the Efimov trimer 
states. For equal scattering lengths, the STM equation 
reduces to that for identical bosons and, in this case, an 
analytic expression for is known 
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167r2Csinh(2?7*) ha^ 
s\v? [£) |a| K,] + sinh^ 77, m 



(for fly = a < 0). (2) 



Here the constants C = 29.62(1) and D = 0.6642(2) have 
been determined from numerical calculations |14l |. The 
analytic expression above exhibits resonant enhancement 
at values a = a,7 = — (e'^/'"')"(DK*)~^ for integer n which 
occur when Efimov trimer states cross the three-atom 
threshold. For 012 7^ 023 7^ ai3, as in our experiment, 
the STM equations must be solved numerically [l^] to 
accurately determine L3. 

Since |ai2|, |a23|, |ai3| > i^dw and ai2,a23,ai3 < 
for fields above 834 G, the data shown in Fig. [2] should 
be well described by the zero-range calculation of L3 in 
Ref. [3l and will allow the first determination of the 
three-body parameters for this three-component ^Li gas 
in the high-field regime. The best fit (solid line) to sub- 
sets of the 30 nK data {B > 907 G) and the 180 nK data 
(B > 975 G) (as described above) is shown in Fig.[21 The 
only free parameters in this fit are and 77*. From the 
fit we obtain = 6.9(2) x lO-^a^^ and r;* = 0.01615^1^?^ 
where the uncertainties - combined statistical and sys- 
tematic - indicate one standard deviation. The zero- 
temperature recombination rate for these three-body pa- 
rameters calculated for fields between 840 and 1550 G is 
shown as the dotted line in Fig. [21 In this model, the 
resonant peak in the recombination rate occurs at 895 G 
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where an Efimov trimer crosses the three-atom scatter- 
ing threshold. The measured three-body recombination 
rate for the 30 nK data also peaks near 895 G, though 
at a significantly smaller magnitude due to the unitarity 
limit. The dashed line shows a "unitarized" recombina- 
tion rate [2§| given by (l/ia + l/issat)^^ to account for 
the observed saturation of the recombination rate in the 
thermal gas to a value L^sat ~ ^Smax/S. 

We can identify the Efimov trimer that crosses thresh- 
old at 895 G with the first excited state of the Efi- 
mov spectrum. The binding energy of this Efimov 
trimer, in the limit of infinite scattering lengths, is Ei = 
^-2-R/sofj2^2^^ _ 55(3) kHz. In this hmit, the ground 
state Efimov trimer has a binding energy a factor of 
g27r/so ^ larger, Eq ^ hx 28(2) MHz. This is the low- 
est state with a binding energy smaller than the van der 
Waals energy scale -Evdw ~ ^^/"^^vdw were the case 
that ai2 — = ^13 ~ a < 0, the first excited trimer 
would cross threshold at a = = —5.0(1) x 10^ and 
the adjacent trimer states would cross threshold at values 
of a either larger or smaller by a factor of e'^/'*" « 22.7. In 
reality, an Efimov trimer of the ^Li system crosses thresh- 
old at 895 G where 012 = — 8584ao, 023 = — 5702ao and 
ai3 = — 2893ao. We identify this trimer as the first ex- 
cited Efimov trimer since falls within this range of 
Uij values. In the high-field region {B > 608 G), the 
ground state Efimov trimer asymptotes to a binding en- 
ergy Eo{a,j = -2140 ao) = hx 24(2) MHz with a width 
7 = (27r) 1.51?-^ MHz (corresponding to a 100^24° ns life- 
time) as calculated from expressions given in Ref. [^. 

One aspect of Efimov physics that can now be tested 
with ultracold atom experiments is the extent to which 
the universal results can be applied when scattering 
lengths are tuned across poles or across nodes. For the 
^Li gas studied here, there are two universal regions 
(122 G < B < 485G and B > 608 G) where the zero- 
range approximation should be applicable, separated by 
a non-universal region where the scattering lengths be- 
come smaller than 2£vdw- We find that the value = 
6.9(2) X 10~'^a^"'^, determined above for the high-field re- 
gion, is close to the value k* = 6.56 x 10~^ Cq ^ previously 
determined from the position of the narrow loss feature 
in the low-field region near 130 G [l3l. Indeed, the re- 
combination rate calculated in the zero-range approxi- 
mation using = 6.9(2) x 10~'^aQ ^ predicts Efimov res- 
onances at 125(3) and 499(2) G in reasonable agreement 
with observations [l|, These resonances occur when 
the ground state Efimov trimer crosses threshold. The 
value of 77* is nearly 1 order of magnitude larger in the 
low field region [lj| than in the high-field regime. The 
significant change in rj^ is likely due to the dramatic dif- 
ference in the binding energy of the deep dimer states for 
the two regions 

Finally, we note that despite the large three-body re- 
combination rate that occurs in this system, we can pro- 
duce quantum degenerate three-component Fermi gases 



in the high-field regime where the three scattering lengths 
arc large, negative and approximately equal. Starting 
from a degenerate two-state mixture of ^Li atoms at 
B ~ 1500 G in trap B, we prepare a three-state mixture 
by applying RF pulses as described above. After equi- 
libration, the three-state mixture has N = 6(2) x 10'* 
atoms per spin state, a temperature T = 50(10) nK, a 
Fermi temperature Tp = 180(20) nK and a degeneracy 
temperature T/Tp = 0.28(6). The difference in mean 
field energies is more than one order of magnitude smaller 
than any other energy scale in the system making it useful 
for future investigations of three-com pon ent Fermi gases 
with SU(3) symmetry [B, [B, 111, 0, llSl ■ 

We arc indebted to E. Braaten and D. Kang for pro- 
viding the computer code to numerically calculate L3. 
We gratefully acknowledge fruitful conversations with E. 
Braaten, B. Esry, J. D'Incao, and D. Kang regarding this 
work. This material is based upon work supported by the 
AFOSR (Award FA9550-08-1-0069), the ARO (Award 
W911NF-06-1-0398) and the NSF (PHY 07-01443). 

During preparation of this manuscript we became 
aware that enhanced recombination in a ^Li ga s near 
900 G was also observed by S. Jochim's group |23l ]. 
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